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We study QCD evolution equations of the first transverse-momentuin-moinent of the naive-time- 
reversal-odd fragmentation functions - the Collins function and the polarizing fragmentation func- 
tion. We find for the Collins function case that the evolution kernel has a diagonal piece same as 
that for the transversity fragmentation function, while for the polarizing fragmentation function 
case this piece is the same as that for the unpolarized fragmentation function. Our results might 
have important implications in the current global analysis of spin asymmetries. 
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^ ■ I. INTRODUCTION 

o . 

Q , The phenomenon of single transverse-spin asymmetry (SSA), An = {(^iSx) — a{—S±))/ {(j{S±) + a{—S±)), defined 
■ as the ratio of the difference and the sum of the cross sections when the transverse spin vector S± is flipped, was first 
observed in the hadronic A production at Fermilab in 1976 as a surprise Since then large SSAs (or other related 
spin effects) have been consistently observed in various experiments at different colhsion energies [2j], such as Sivers 
and Collins asymmetries in the semi-inclusive hadron production in deep inelastic scattering (SIDIS) and in hadronic 
collisions 0, |4| , as well as the large cos(20) anomalous azimuthal asymmetry in back-to-back dihadron production in 
' e^e^ annihilation Q. 

To understand all these non-trivial and interesting spin effects and to explore the physics behind these asymme- 
O ■ tries, two QCD-based approaches have been proposed and widely applied in the phenomenological studies [6l-[l^: 
D I the collinear twist-three factorization approach [Ij] and the transverse momentum dependent (TMD) factorization 
i-C ■ approach [l3 - [l8j . These two approaches apply in different kinematic domain, and have been shown to be equivalent 
in the overlap region where they both apply; thus provide a unified QCD description for these spin effects 19]. In the 
. collinear twist-three factorization approach, the spin effect depends on certain twist-three multi-parton correlation 
I> functions. On the other hand, in the TMD factorization approach, the spin effect could be described in terms of TMD 
On , distributions and fragmentation functions. These twist-three multi-parton correlation functions are closely related 
to the TMD functions. For example, so-called Efremov-Teryaev-Qiu-Sterman quark-gluon correlation function is the 
first transverse-momentum-moment of the Sivers function [18j . 

The spin effect could be generated from either the spin correlation in the parton distribution functions, among which 
' the Sivers [l^ and Boer-Mulders 17] functions are the important examples; or the spin correlation in the fragmentation 
functions, among which the Collins function [isj and polarizing fragmentation function [l6j are the important ones. 
Although all these four functions are naive-time- reversal-odd (T-odd), they have very different universality properties. 
For both Sivers and Boer-Mulders functions, it has been shown that they differ by a sign for the SIDIS and Drell- 
Yan (DY) processes [20|. On the other hand, both Collins function and the polarizing fragmentation function have 
been argued to be universal between different processes (2l| - [26| . The different universality properties are connected 
with the non-trivial initial- and final-state interactions between the active parton and the target remnant p7| , whose 
d ' interesting consequences remain to be tested in the future experiments 0, [26| . 

On the fragmentation side, both Collins function and polarizing fragmentation function have been widely used in 
describing the spin asymmetries observed in the experiments. Collins function describes the transversely polarized 
quark jet fragmenting into an unpolarized hadron, whose transverse momentum relative to the jet axis correlates with 
the transverse polarization of the fragmenation quarks. It has been believed to be responsible for the azimuthal asym- 
metry observed in SIDIS 0, and the cos(2(/)) asymmetry observed in back-to-back dihadron production in 
e~^e~ annihilations On the other hand, polarizing fragmentation function describes the distribution of transversely 
polarized hadron in an unpolarized quark, through a correlation between their relative transverse momentum and the 
hadron transvers e sp in vector, which have been believed to be responsible for the Hyperon polarization observed in 
the experiments |28j . 

The first transverse-momentum-moment of the Sivers and Boer-Mulders functions corresponds to twist-three quark- 
gluon correlation functions Tq^p{x,x) and Tj'^(x,x) fi^, for which the QCD evolution equations have been studied 
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by various authors [29|-|32|. On the other hand, the corresponding fragmentation correlation functions connected to 
the first transverse-momentum-moment of the Cohins function and the polarizing fragmentation function have been 
identified only recently [IE HB] • The purpose of our paper is to derive the QCD evolution equations for these relevant 
fragmentation correlation functions. QCD evolution equations are important in the sense that they control the energy 
dependence of the associated spin observables, also they enable us to evaluate the higher-order corrections to the 
spin-dependent cross sections systematically. 



II. EVOLUTION EQUATIONS 



To start, we recall that the eight leading-twist TMD fragmentation functions could be defined from the following 
correlator 



Zh J (2Try 



^{0\C^i;mPSX){PSX\m^l\0) 



(1) 



where a factor l/27Vc from the average over the spin and the color of the fragmenting quark is suppressed. P is the 
momentum of the final-state hadron with spin S, which has a transverse component p± relative to the momentum k 
of the fragmenting quark. We choose the hadron to move along the +z direction, and define two light-like vectors: 



[l+,0-,Oj 



= [o+,i-,oJ 



(2) 



The momentum fraction Zh = P'^ /k'^ , and k±_ = —p±/zh- To ensure gauge invariance, we have explicitly written out 
the gauge link — V exp [~ig Jg°° dXn ■ A{£^ + Art)) . 

The correlator A{zh,p±^) could be expanded as follows [iIgI. [TtI 



Mzh,p±) 



Di{zh,pi)r^ + XhGiL{zh,piWT^ + Hiizh,p\)iaap-f^n"S'l 



+GiT{zh,Pi_)—rj — 7^ + >^hH^L{zh,p^) 



Ml 



-l(Tal3l n 
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where A/i is the hadron helicity and Alh is the hadron mass. Out of the above eight TMD fragmentation functions, 
Di{zhTp\), GiL{zh,p'±), and Hi{zh,p']J arep_L-even functions and correspond to the unpolarized, longitudinally and 
transversely polarized distributions in the fragmentation. They are related to the leading-twist collinear fragmentation 
functions after integrating over p±, for which the QCD evolution equations have been well-known [33, 34]. All other 
five TMD fragmentation functions are pj_-odd functions in the sense that they vanish if integrate over p± directly. 
However, for D^rp{zh,p'±), H^{zh,p'±), GiT{zh,p'±), and H^j^lzhjp"^), their integral over p± after first weighted by 
p\ (called "first p^-moment") lead to the twist-three collinear fragmentation correlations; for H^{zh,p']_), the non- 
vanishing integral needs to be weighted by even higher power p± and actually corresponds to twist-four fragmentation 
correlation. 

In this paper, we are particularly interested in deriving the QCD evolution equations for the first transverse- 
momentum-moment of Hi{zh^p\) and Dj^rp{zh^p\): H^{zh,p\) is the Collins function, and D^rp[zii,p']_) is the 
polarizing fragmentation function, which have been the main focus in the phenomenological studies. Their first 
transverse-momentum-moments have been identified recently, and given by the following fragmentation correlation 
functions fHIH 



{Tr 7^(01 



f{z) = 



27r 



nnS± a 



^{n\PSX){PSX\m\0) + h.c], (4) 
i>ir)\PSX){PSX\m\0) + h.c.}, (5) 
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where D" = 9" — igA'J_ is the covariant derivative, F"^ is the gluon field strength tensor. H{z) and T{z) are related 
to the Collins function H^{zh^p'\) and polarizing fragmentation function D^j,{z)^,p'\) as follows: 



H{z) 
f{z) 



d P^^rrHi {z,P±), 
d'p^^^Dip{z,pl). 



(6) 
(7) 



According to Refs. [25|, |26|, |35| , the above defined one-variable fragmentation correlations belong to the more general 
twist-three two-variable fragmentation correlations, which are defined as 



where fc+ 



Hd{zi,Z2) 
Td(zi,Z2) 

k1 



Z1Z2 f dr] 
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Z1Z2 



2tt 2tt 



2-K 27r 



-e^'=.>"-{Tra"^n^(0|ii^i„(?7-)7/^(r)l^^^)(^^^^lV^(0)|0) + /i.c.}, (8) 
^-K^- 1 |Tr me^^^'^'^iD^^ {rj^ ) ^(r ) 1^^^^) {PSX\^{0)\0) + h.c. } , (9) 



fc^ with k'l 



I z\ and fc^ ~ P+/z2- Similarly, one can define the corresponding _F-type 

fragmentation correlations Hf{zi,Z2) and Tf{zi,Z2) by replacing Z?" by qF'^ °' in Eqs. ([8|) and ([9|). By using the 
equation of motion, _D-type and i^-type functions are related to each other [Tsl , [25l . l36j 



Hd{zi,Z2) = PV 

f,,(zi,Z2) = PV 




Hpizi, Z2) + S 



1 
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^2 / Z2 
Zl, 



H{Z2), 



TF{z^,Z2)+&[---]-nZ2), 
Zl Z2 Z2 



(10) 
(11) 



where PV stands for principle value. The above equations show that Hd and Tu are more singular than their 
corresponding _F-type functions Hp and Tp at z\ = Z2. Moreover, recent study has shown that Hp{zi, Z2) and 

Tp(zi,Z2) vanish when zi — Z2 [2^ [23|. Therefore, it is more convenient to use the set of ^H(z), Hp(zi, Z2)^ and 
|r(z), Tp(zi, 2:2)! when we derive the evolution equations, as we will follow below. 






(a) 

FIG. 1: Feynman diagram representation: (a) the first transverse-momentum-moment of TMD fragmentation functions, (b) 
evolution contribution from itself: I 



correlations: (.q = P/z, iq-^ = P/zi, and £g 



P/ z -\- £q^ , (c) evolution contribution from the two-variable _F-type fragmentation 



The fragmentation correlations H{z) and T(z) could be represented by the same cut forward scattering diagrams 
as sketched in Fig. [ija), but with different cut vertices which are used to connect the operator definition of the 
fragmentation correlations to the cut forward scattering Feynman diagrams. The standard way to derive the cut 
vertex is to express the operator definition of the correlation functions in terms of hadronic matrix elements of quark 
and gluon operators in momentum space, for details, see [2^. The cut vertices for H{z) and T{z) as represented in 
Fig. [ija) are given by 



k^j-n, (12) 

^nnS^k^^ (13) 







P+ 

z 






p+ 

z 



where k± is the transverse momentum of the fragmenting quark with respect to the final hadron momentum. 
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Let's now explain how to derive the evolution equations. We will work in the light-cone gauge n- A = 0. In order to 
derive the evolution equations and evolution kernels from the operator definitions of the fragmentation correlations, 
one needs to compute the perturbative modification to these functions caused by the quark-gluon interactions in QCD. 
For both H(z) and T(z), the perturbative modification could come from either these correlation functions themselves, 
as shown in Fig.[ljb); or the corresponding F-type fragmentation correlations, as shown in Fig.[IJc). Since both H{z) 
and T{z) correspond to the operator ~ {0\d±ip\PSX){PSX\'ijj\0), i.e., the partial derivative in the quark field, in 
order to calculate the contribution from themselves as shown in Fig.[ljb), one has to perform coUinear expansion. In 
other words, one should assume £q ~ P/ z + , and the linear in iq^ expansion term when combined with the quark 
field from the top blob will lead to the fragmentation correlations H{z) and T(z). On the other hand, to calculate 
the contribution from the -F-type fragmentation correlations as shown in Fig. [Ijc), one has to insert a gluon in 
the Feynman diagram and then convert A" to field strength through a partial integration, which leads to the 
-F-type fragmentation correlations. In the calculation of A±^ contributions, since no collinear expansion is involved. 



one could set all the parton momenta as collinear to the hadron: Iq 
up, the perturbative modifications could be written as 



P/z, = P/zi, and £g 



To sum 



dH{zh,n'^) 



df(zh,fJ.'^) 



dz 



_1 d_ 



Tr [i • P/C(fc, £q « P/z + Iq^ )]^ /i^) + / dzdziPV - 




Tr \i 7"7 • PK^(k, - P/^', ^,1 = P/^i)] Hf{z, zi, a^'), 



(14) 



^TT^Tr [7 • Pe""^^"/C„(fc,€, = P/zJq, - P/zi)] fpiz.zi.i?), 



(15) 



where /i is the factorization scale, IC{k,iq Pj z \- £q^) and ICa{k,£q — P/z,lq^ = P/zi) are the hard partonic part 
calculated from Fig.[TJb) and (c) without the top blob, respectively. Note we use the same symbols /C(fc, £q k, Pj z+£q^ ) 
and ICa{k,£q = P/z,£q-^ = P/zi) for both H{zh,n^) and T(z;i,/i^) for simplicity, they are different in the calculations 
as shown below. 





(c) 

FIG. 2: Contribution from the first transverse-momentum-momcnt of the TMD fragmentation functions themselves: (a) real 
contribution, (b) and (c) virtual contributions. 



To the leading order in strong coupling constant as, the contribution from the fragmentation correlation functions 
themselves as in Fig.[ljb) are given by the Feynman diagrams in Fig.O [2][a) is the real contribution, while [2][b) and 
(c) are the virtual contributions. These diagrams are the same as those when one calculates the evolution kernel for 
the leading- twist unpolarized collinear fragmentation function [H,!!^], but the actual calculations are very different. 
As we have explained above, in our current calculations, a collinear expansion is needed to pick up the linear in £q^ 
terms which lead to H{z) and T{z)] on the other hand, there is no collinear expansion involved in the calculations 
for the leading-twist collinear functions. Following the collinear expansion as specified in the first terms of Eqs. (fT4|) 
and (|15|) . we obtain the contributions from the real diagram Fig. [2l^a): 

where z = Zh/z. For virtual diagrams Fig. [DJb) and (c), since now £q — k, thus the £q^ expansion is in fact an 
expansion of k±. On the other hand, the cut vertices for both H{z) and T{z) in Eqs. and depend linearly 



5 



on k±, thus after a direct expansion over k± in the cut vertices, one could set ah k±_ = £q^ 
resuhs are 



b) + (c 



dTizh,n )lFig.[lb)+(c) = 
Combining above real and virtual contributions, we obtain 



di?(z;„/x2)|p.g [2] = /" ^^Cf 



dk^ as 



k\ 27r 



Cf 



^ dz'l + Z'2 , 



afterwards. The final 

(18) 
(19) 



df{zh,^l^)\^.^^ 



fc2_ 27r 



z 

dz 

z 



22 3„_, 

(l-z)+ 2 



f{z,^?). 



(20) 
(21) 



Let's now consider the contribution from the i^-type fragmentation correlation functions. To the leading order 
in Q!s, the relevant Feynman diagrams are shown in Figs. [3] and HI the real contributions in Fig. [3] and the virtual 
contributions in Fig. U) As we have mentioned above, calculating the A± contribution (finally related to F-type 



(a) 

FIG. 3: Contribution from the _F-type fragmentation correlation functions: real diagrams. The "mirror" diagrams for which 
the additional gluon attaches on the left of the cut are not shown, but are included in the calculations. 

fragmentation correlation functions) does not involve collinear expansion in the light-cone gauge, we could set all the 
parton momenta collinear to the final hadron: 






P/z, 



^91 



~ ^qi • 



(22) 



The calculations following the formalism in the second terms of Eqs. and are also straightforward. For the 
real diagram contributions, we collect the terms through the color factors: Fig. [Sja), (b), and (c) have color factors 
Cf, Cf — CaI"^, and Ca/2, respectively. The final results are 



k^ 27r 



dz f dzi 



— / — PV B{zh,z,zi)HF{z,zi,iJ, ), 



dk'^ as f dz f dzi 



kl 2n 

where the kernel B(zh, z, zi) and B'{zh, z, zi) are given by 



— / — PV -J B'izh,z,zi)TF{z,zi,n ), 



B{zh,z,zi) = Cf 
B'{zh,z,zi) = Cf 



2zh 

z 



1 



Zh_ _ Zf^ 

Z\ z 



Ca 
2 



2zfe zliz^ + z\) - zzi{z + zi) 
z {zi - z){zi - Zh)z 



Zl Zl Z ZZ\ 



Ca 
2 



{zzh + ziZh- 2zzi){zzi + zl) 



[zi - z){zi - Zh)z'^ 



(23) 
(24) 

(25) 
(26) 



Finally let's consider the virtual contributions from the _F-type fragmentation correlation functions as shown in 
Fig. m It is important to realize that for all the diagrams (a)-(e) in Fig. HI we have (follow Eq. ([221) ') 



k = £q = P/z, 



(27) 



which has no transverse component, i.e., fc_L = 0. Note the cut vertices used to define both H{z) and T{z) depend 
linearly on k±, see Eqs. (IT^ and p3p . Thus when k± = 0, they vanish. In other words, all these virtual diagrams 
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do not contribute. Thus the perturbative modifications for H[zh,^i^) and T{zh,^i^) receive contributions from only 
Figs. [5] and [21 Adding them up, we obtain 



dfc^ as f dz 



2n dkl as f dz 



A{z)H{z,n'')+ I ^FY 1^^) B{zh,z,Zl)HFiz,Zl,^i^) 



A'iz)f{z,fj,^)+ I ^pv|^j^j S'(z^,z,zi)fj.(z,zi,^2) 



(28) 
(29) 



Differentiate both sides of above equations with respect to In/x^, we obtain the scale evolution equations for H{zh, /i^) 









(a) 



(b) 



(c) 



(d) 



(e) 



FIG. 4: Contribution from the _F-type fragmentation correlation functions: virtual diagrams. The "mirror" diagrams for which 
the additional gluon attaches on the left of the cut are not shown, but are included in the calculations. 



and T{zh, p?) as 



dH{zu,ii^) 

df{z,,,fi^) 
9 In ^2 




B{zh,z,Zi)Hp{z,Zi,n'^) 
B'{zh,z,zi)fF{z,zi,n^) 



where B{zfi, z, zi) and B'{zh, z, zi) are given in Eqs. ([25)1 and (|26l) . and A{z) and A'{z) have the following forms 



A{z) = Cf 
A'{z) = Cf 



2z 



(l-z) 
1 + z' 



+ ^<5(l-i) 



(30) 
(31) 

(32) 
(33) 



.(l-z) 

Eqs. (pO| and ((3T|) are the main results of our paper. A few comments about these results are provided: 

• The evolution equations derived here for both H{zh,n^) and T{zh,n'^) are not a close set of equations, as 
stand in Eqs. (pO|) and ([3T|) : the evolutions depend on the diagonal pieces H{z, /i^) and T{z, /i^), as well as the 
off-diagonal pieces Hf{z, zi, p?) and Tf{z, zi, y?). 

• It is interesting to notice that the evolution kernel A{z) is the same as that for the transversity fra gme ntation 
fimction [34], while the kernel A'{z) is the same as that for the unpolarized fragmentation function [33| . 

• As shown in Refs. [23l|. both Hf{z, zi, y?') and TV(^,^i,M^) vanish at z = zi. This might imply that the 
off-diagonal pieces could be small, thus the evolution of H{zh, fJ.^) might be close to that of transversity, while 
the evolution of T{zh, /i^) might be close to that of unpolarized fragmentation function. If this were true, it will 
have important consequences on the current global analysis of the spin asymmetries @, Q • Of course, whether 
the off-diagonal terms play a less important role in determining the evolution of the diagonal terms needs to be 
tested from experimental data through global analysis, such as those done in [ol. [Toj. 

• It is also important to realize that there is no gluon Collins function, thus H{zh, n^) does not receive contribution 
from gluon part. On the other hand, there could be gluon polarizing fragmentation functions (ssj . from which 
the corresponding gluon fragmentation correlation functions could be defined. Thus there could be contributions 
from the gluon part to the evolution of T{zh, /i^), and these contributions are not studied here. 
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III. CONCLUSIONS 



We have derived the QCD evolution equations for the first transverse-momentum-moment of the naive-tinie-reversal- 
odd transverse momentum dependent fragmentation functions: the ColHns function H^{zfi,p']_) and the polarizing 
fragmentation function D^rp{zh,p]_)- These first transverse- momentum- moments correspond to twist-three fragmen- 
tation correlation functions denoted as H{zh,fi'^) and T{zh,fJ.'^)- We calculate in light-cone gauge the order of 
evolution kernel for the scale dependence of both H{zh,fJ.'^) and T{zh,iJ.'^)- We find that the evolution of both 
fragmentation correlation functions receives contributions from themselves, as well as from the F-type two-variable 
fragmentation correlation functions Hf{z, zi, ^'^) and Tf(z, zi, p?). We find for H{zh,^^) that the diagonal piece 
in the evolution kernel is the same as that for the transversity fragmentation function, while for T(zft,/i^) that the 
diagonal piece is the same as that for the unpolarized fragmentation function. Since the off-diagonal pieces involve 
the -F-type fragmentation correlation functions which vanish at z = zi , thus they might play a less important role. If 
this were true, it will provide important consequences in the current global analysis of spin asymmetries. 
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